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We give a general method to find exact cosmological solutions for scalar-field dark energy in the 
presence of perfect fluids. We use the existence of invariant transformations for the Wheeler De Witt 
(WdW) equation. We show that the existence of a point transformation under which the WdW 
equation is invariant is equivalent to the existence of conservation laws for the field equations, which 
indicates the existence of analytical solutions. We extend previous work by providing exact solutions 
for the Hubble parameter and the effective dark-energy equation of state parameter for cosmologies 
containing a combination of perfect fluid and a scalar field whose self-interaction potential is a power 
of hyperbolic functions. We find solutions explicity when the perfect fluid is radiation or cold dark 
matter and determine the effects of non-zero spatial curvature. Using the Planck 2015 data, we 
determine the evolution of the effective equation of state of the dark energy. Finally, we study the 
global dynamics using dimensionless variables. We find that if the current cosmological model is 
Liouville integrable (admits conservation laws) then there is a unique stable point which describes 
the de-Sitter phase of the universe. 
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1. INTRODUCTION 

The discovery of the accelerated expansion of the universe (see jjJ and references therein) has opened a new window 
in cosmological studies. Indeed, the underlying physical process responsible for this phenomenon is considered as one 
of the fundamental problems in cosmology. Within the framework of general relativity, scalar fields provide possible 
dark energy models which can describe, but not so far explain, this acceleration. Scalar field models require the choice 
of a self-interaction potential V((f>) for the scalar field <j>. The considerations of a specific potential V{<j>) is done by 
an adhoc requirement (ansatz). In this manner various candidates have been proposed in the literature, such as an 
inverse power law, exponential, hyperbolic and the list goes on (for review see 0 and references therein). One such 
ansatz, which we shall consider in the following, has been done by Rubano and Barrow Q (see also i) who found 
that if the scalar field behaves as a perfect fluid then the potential V((j>) oc sinh p (g^), where the constants q,p are 
given in terms of observable cosmological parameters, namely the dark-energy equation of state (EoS) parameter and 
flmO • 

Obviously it is important to proposed potentials which are realistic and at the same time lead to exactly soluble 
models in the FLRW spacetime. The reason for employing the Noether symmetries is that Noether symmetries 
provide us with the Noether integrals which facilitate the analytic solution of the field equations. As an example, in 
extended theories of gravity, where the Birkhoff theorem is not guaranteed, the Noether approach provides a means of 
describing the global dynamics using the first integrals of motion [5], Furthermore, besides the technical possibility of 
reducing the dynamical system, the first integrals of motion always give rise to conserved currents that are not only 
present in physical space-time but also in configuration spaces (see discussion in 0 )- In space-time such currents are 
linear (momentum, angular momentum etc.) but in configuration space the conserved quantities appear as relations 
between dynamical variables The latter implies that in the configuration space (minisuperspace), the first integrals 
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are considered as ’’selection rules” for scalar field potentials and coupling functions in the case of modified gravity 
theories. The above features have inspired many authors to propose the admittance of a Noether (point or dynamical) 
symmetry by the field equations as a selection rule for dark-energy models, including those of modified gravity [Ml- 

In particular^ it has been shown that the dynamical Noether symmetries are associated with the Killing tensors 
of the space [li Hi and that the Noether point symmetries are related to the homothetic algebra of the space (see 
[22i ] and references therein). Therefore the requirement of a Noether symmetry is indeed a geometric demand hence 
independent of the particular dynamics. 

In scalar field cosmology the system develops in minisuperspace whose geometry is defined by the field equations. 
Therefore the requirement that the field equations admit a Noether symmetry becomes, according to the above results, 
a geometric requirement on a geometry which is inherent to the system. 

In this work by using a more general geometric criterion, i.e. by employing the Lie point symmetries of the Wheeler- 
DeWitt (WdW) equation, we extend the work of Rubano and Barrow Q to a general family of hyperbolic scalar-field 
potentials V(<p). 

In section [21 we give the basic theory of the scalar-field cosmology in a FLRW spacetime. The basic definitions and 
results from the Lie and Noether point symmetries of partial differential equations and the application in the WdW 
equation are presented in section [3] In section [IJ we consider our cosmological model, which includes a hyperbolic 
family of scalar field potentials with a perfect fluid with constant equation of state parameter w m . We study the 
existence of Lie point symmetries of the WdW equation, where we find that for our model the WdW equation admits 
Lie point symmetries if the free parameters of the potential and the parameter w m are related. In section [5] we apply 
the Lie point symmetries of the WdW equation in order to construct invariant solutions of the WdW equation and 
exact solutions of the field equations. Also, in section [5] we perform a dynamical analysis by studying the fixed points 
of the field equations in the dimensionless variables for the general model and we show that when the cosmological 
model is Liouville integrable (the model admits conservation laws), there is a unique stable point which describes the 
de Sitter universe. Finally, in section [7] we draw some conclusions. 


2. SCALAR-FIELD COSMOLOGY 


We start with the Friedmann-Lemaitre-Robertson-Walker (FLRW) spacetime with line element (c = 1) 

1 


ds 2 = —dt 2 + a 2 (t ) 


(1 + f X 2)2 


( dx 2 + dy 2 + dz 2 ). 


(1) 


The total action of the field equations is written as 


S = Seh + S# + S m , (2) 

where Seh = / dx 4 ^/—gR is the Einstein-Hilbert action, R is the Ricci scalar of the underlying space, S0 is the action 
of the scalar field 


S</, = / dx^yf^g 


+ v^) 


( 3 ) 


and S m = J dx^yJ—gL m is the matter term. We assume that 0 inherits the symmetries of the metric ([!]) therefore 
0(t) and consequently 0 ; „ = 05° where 0 = 

From the action ([21), we have the Einstein field equations Hi 


P/ll/ 


( 4 ) 


where k = 8nG = 1, R is the Ricci tensor and is the total energy momentum tensor given by +1)^(0). 

T^ v is the energy-momentum tensor of baryonic matter and radiation and T Ml/ (0) is the energy-momentum tensor 
associated with the scalar field 0. Modeling the expanding universe as a fluid (which includes radiation, matter and 
DE) with 4—velocity proper isotropic density p m and proper isotropic pressure P m gives = — P g+ (p + 
P)u^u„, where p = p m + p ^ and P = P m + P^. The variable denotes the energy density of the scalar field and P^ 
is the corresponding isotropic pressure. Moreover the parameters {p < / > ,P ( j > ) of the scalar field are given by 

p 0 = ^0 2 + K(0) , P 0 = i0 2 - V(0). 


( 5 ) 





3 


For the FLRW spacetime, for comoving observers (tt M = S(f), the Einstein field equations J4]) are 


H 2 


3 {pri 


K 

712 


and 


3 H 2 + 2 H = - K {P m + P 0 ) - —, 

a z 


( 6 ) 

(7) 


where H(t) = a/a is the Hubble function. 

Furthermore, assuming that the scalar field and matter do not interact, we have the two conservation laws 


Pm + 3H(p m + P m ) — 0 


( 8 ) 


p</> + 3H(p,j J + P^) — 0 


(9) 


while the corresponding equation of state (EoS) parameters are given by w m = P m /p m and w$ = P ^/In what 
follows we assume a constant w m , so that p m = p m oa~ 3< ' 1+Wm ^ ( w m = 0 for cold matter and w m = 1/3 for radiation), 
where p m o is the matter density at the present time. Generically, some high-e nerg y field theories suggest that the 
dark energy EoS parameter may be a function of cosmic time (see, for instance, 

Replacing (|5j in © we have the Klein-Gordon equation 


2 |). 


+ 3 Hct> + = 0 

where Furthermore the corresponding dark energy EoS parameter is 


( 10 ) 


P* (/ ) 2 / 2 ) - VW 

W * P4> {j> 2 / 2) + V(0) 


(ID 


which means that when then ft < V{<f>). On the other hand, if the kinetic term of the scalar field is 

12 

negligible with respect to the potential energy, i.e. <C V (</>), then the equation of state parameter is w</, ~ — 1. 

From the above analysis it becomes clear that the unknown quantities of the problem are a(t), <j>(t) and V(<j>) 
whereas we have only two independent differential equations available namely Eqs. © and m- Therefore, in 
order to solve the system of differential equations we need to make an additional assumption (ansatz). This usually 
concerns the functional form of the scalar field potential V (</). In the literature, due to the unknown nature of DE, 
there are many forms of this potential which describe differently the physical features of the scalar field (for instance 
see Hi, SMI)- 

As far as the exact solution of the field equations ©, ® and m is concerned there are few solutions with spatial 
curvature [33|, Hi] and even fewer solutions are known for a perfect fluid and a scalar field jUl, l35l - [38[ . A special 
solution for a spatially fiat FLRW spacetime ( K = 0) which contains a perfect fluid with a constant equation of state 
parameter P m = (7 — 1) p m and a scalar field with a constant equation of state parameter = 7 ^ — 1 = P^/p^, has 
been found in Q. Specifically in Q it has been shown that under these assumptions one solves the field equations 
and finds the potential V ( (f >): 


V ((/>) = 3H 2 (1 - O m0 ) (l - 


f 1 ^mO 

\ ^mO 


\ ~ 2± - 
\ 7 — 


sinh 


( ,/52 7# ( ^' w ) 


2 7 (f, 
7-7 $ 


( 12 ) 


Evidently, this solution is a special solution, in the sense that it exists for specific initial conditions, for example with 
Wcf, (z) constant. 

In the following we consider a spatially-flat FLRW spacetime with a perfect fluid P m = (7 — 1) p m and a scalar 
field and assume that the potential has the generic form 


V (</>) = Vo [a cosh (p<j)) + /3 sinh (p </>)] 9 , (13) 

where Vo,a,/3,p,and q are constants. This potential (fT51) is a generalization of (fl^l) and our aim is to determine 
exact solutions of the field equations for a particular relation between the constants p, q of the potential (1131) and the 
barotropic parameter 7 of the perfect fluid. 
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3. PRELIMINARIES 

In this section, we show that if the WdW equation admits Lie symmetries which form an Abelian Lie algebra, then 
the WdW equation is Liouville integrable; that is, the field equations can be solved by quadratures. 


3.1. Lie and Noether point symmetries 


A Lie symmetry of a differential equation H = H(x l , u A , , u^) is the generator of the one parameter point 

transformation which leaves invariant the differential equation H. That means that if X =C(x k , u B )dj+ri A (x k , u B )d rj A 
is a Lie symmetry for H, then there exists a function A such that the following condition holds jihil . [40| 

X [ 2 1(U) = A H , mod H = 0 (14) 

where X.^ = X + r]fd u A + pfjd u A , is the second prolongation vector 1 of X. 

The importance of Lie symmetries is that each symmetry can be used to reduce the number of dependent variables. 
Solutions which follow from the application of Lie symmetries are called invariant solutions. 

For differential equations which arise from a variational principle there exists a special class of Lie symmetries, the 
Noether symmetries. Noether symmetries are Lie symmetries which leave the action integral invariant. According to 
Noether’s theorem to each Noether symmetry there corresponds a conserved Noether current. 

The condition for a Noether symmetry is that there exists a vector field A 1 = A l (x z ,u ) such that the following 
condition is satisfied: 


X [1] L + LDiC = D t A l . 


(15) 


The corresponding Noether current I 1 is defined by the expression 


r = i k 



a Ai 

’ s ? + ' 4 ' 


( 16 ) 


and it is conserved, that is, it satisfies the relation D,P = 0 (39| . 

As we discussed above, the method of using the Noether symmetries of the cosmological field equations has been 
applied by many authors in scalar-field cosmology, in f(R) gravity and other modified gravity theories. Recently 
|23l l41j , it has been proposed that the cosmological model will be determined by the existence of Lie symmetries of 
the WdW equation of quantum cosmology. 

This selection rule is more general than that imposed by the Noether symmetries of the field equations, because, as it 
has been shown in [4lj| . the WdW equation is possible to admit Lie point symmetries while the Lagrangian of the field 
equations does not admit Noether point symmetries. In the following we discuss the application of Lie symmetries 
in the WdW equation. Specifically, we discuss the reduction process and we show how to construct Noetherian 
conservation laws for a conformally related Lagrangian of the field equations from the Lie point symmetries of the 
WdW equation. 


3.2. Minisuperspace and invariant solutions of the WdW equation 

The Lagrangian of the field equations in minimally coupled scalar field cosmology in a spatially flat FLRW spacetime 
with a perfect fluid with a constant equation of state parameter P m = (7 — 1) p m is 

L (a, a, </>, cfj = —3aa 2 + — a 3 (j) 2 — a 3 V (ft) — p m oa _3 ^ 7_1) - (17) 

The field equations are the Euler-Lagrange equations of m with respect to the variables (a, 4 >) and are equations 
m and GOD. As the Lagrangian is independent of time, we also have the Hamiltonian constraint (0), which, in terms 
of the momenta p a = §§, P 4 , = becomes 

- -^pI + 7 ^ 3 Pi + a 3 V ((/>) + pm 0 a~ 3{l ~ 1) = 0 . ( 18 ) 


1 Where 


= Da 


) 


— „• uDn (£ k ) and Do is the total derivative. 

i k 3 V ' S / 1 
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Finally, the field equations are equivalent to the following Hamiltonian system: 


a = - 7 -Pa , 
o a 


1 s 

= —P0 , = -O V,0, 


Pa = - - 3a 2 F (0) + (3y - 3) pmoa 37+2 = 0. 

The WdW equation is the Klein Gordon equation which is defined by the conformal Laplacian operator. The 
general conformal Klein Gordon equation is: 


AT 


n — 2 
4 (n — 1) 


R(x k )^ + V eff (z fc )T = 0, 


(19) 


where A = -j= (vlfffgfj) is the Laplacian operator, is the metric of the space and n = dim cjij. We brake the 
Lagrangian m in two parts. The kinematic part which we consider as a Riemannian space with dimension n = 2, 


and line element 


ds 2 = —6ada 2 + a 3 d(p 2 (20) 

which we call the minisuperspace, and the dynamic part which is defined by the potential V e f / (a, <p) = 

2a 3 [V (</>) + Pmoa -37 ] • Specifically, the minisuperspace is a 2-dimensional Lorentz manifold whose coordinates are 
the scale factor and the scalar field. The symmetries of this space are related to the Lie and the Noether symmetries 
of the dynamical field equations. The metric of the minisuperspace is defined by the kinematic part of the Lagrangian 
(fT71) for the dynamics of the field. The rest of the Lagrangian is considered to be the "effective potential” of the 
dynamical system that is defined by the gravitational and scalar fields. Therefore, using the minisuperspace (1201) the 
WdW equations becomes 


AT + 2a 3 [V (</>) + p m0 a" 37 ] * = 0, 


( 21 ) 


where the Laplacian operator A is given by 

A=-—(— — 

6 a \9o 2 da 

In (23[, it was proved that the Lie point symmetries of equation (flTH) are related to the conformal algebra of the 
minisuperspace gtj. More specifically, it has been shown that: 

A. The general form of the Lie point symmetry vector is 


1 d 2 
a 3 d(f> 2 


( 22 ) 


X = f ( x k ) d, + 


(2 ~ n) 
2 


ij) T + aoT 


d'S', 


(23) 


where £* ( x k ) is a conformal Killing vector of the minisuperspace, with conformal factor ip (x fc ) . 

B. The Lie point symmetry condition which constrains the potential is L^V e ff + 2ipV e ff = 0. 

We require now that equation m admits as Lie point symmetry the vector (1231) . Then under the coordinate 
transformation x l y l so that (a; fc ) di —)• dj , the Lie symmetry vector (E3l) becomes 


x = dj + 


2 — n 
2 


ip T + a 0 T 




(24) 


There exist two equivalent methods to reduce the WdW equation by means of the symmetry vector (l24l) . 
a) In the first method we calculate the zero-order invariants from the Lagrange system (6 ^ J), 

dy b dy J dT 

— = ~T = (^ip + a 0 ) T’ 

which turn out to be 


y b , T (y b ,y J ) = T (y b ) exp 



-ip + a 0 



(25) 


(26) 
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b) The second method is to write the Lie point symmetry as a Lie Backlund symmetry. The Lie point symmetry 
PI is equivalent to the contact symmetry 


X= 

from which we obtain the differential equation 


2 — n 


ip + ao ) 'F ) dq, 


2 — n 


*F j - 2 ‘" tp + ao J *F = aid'. 

We set ao + ai = Q o and find that the solution of the reduced equation is 

'2 — n 


^ (y b ,y J ) = $ (y b ) exp J [-J-+ Qo J dy 


(27) 


(28) 


(29) 


from which it follows again that the coordinate y J is factored out from the solution of the wavefunction ’I' (y fc ,j/‘ 7 ) . 

In the WKB approximation, *F ( x k ) ~ e lS ( x ) the WdW equation reduces to a (null) Hamilton-Jacobi equation. 
The latter can be seen as the Hamilton-Jacobi equation of a Hamiltonian system moving in the same geometry under 
the conformal Laplace operator of the WdW equation and with the same potential. Specifically, the WdW equation 
m in scalar field cosmology provides the null Hamilton-Jacobi equation: 


12a \da J 


1 

2^3 



+ a 3 V (</>) + p m0 a 3(7 15 = 0. 


(30) 


Furthermore, in [23| it was also shown that the symmetries of the WdW equation can be used in order to find 
Noether point symmetries for classical particles. However, the null Hamilton-Jacobi equation is separable if the n- 
dimensional Hamiltonian system admits n conservation laws <Jq (symmetries) i.e. n corresponding Noether symmetries 
which are independent and in involution, i.e. {<F/, d>^} = 0 where .} denotes the Poisson bracket. If this is the 
case, then the Hamiltonian system is Liouville integrable [42j| . That means that it is possible for the WdW equation 
to admit an invariant solution and at the same time the classical Hamiltonian system to be not integrable. Therefore, 
in order for the WdW equation to admit an invariant solution and the Hamiltonian system to be Liouville integrable, 
the n- dimensional WdW equation must admit at least n — 1 independent Lie point symmetries, A/, which form an 
Abelian Lie algebra. If this is the case, the zero-order invariants of these n — 1 Lie point symmetries will give the 
solution of the WdW equation in the form 


'F ( x n ,x J ) = <F (x n ) exp 


‘n—1 

E 

,j= i • 


2 — n 


ip — Qj) dx 


(31) 


where Qj are constants, J = 1, 2,..., n — 1, and the function $ ( x n ) satisfies a linear second-order ODE. That is, when 
the field equations are Liouville integrable by Noether point symmetries then there exists a coordinate system where 
the WdW equation admits n oscillatory terms in the solution and vice versa. It is important to note that this result 
is more general and includes the one given in (13] when ip ( 'x k ) = 0; that is, if one considers the Killing algebra of the 
minisuperspace only. We conclude that for the reduction/solution of the WdW we may consider directly the Lie point 
symmetries of the WdW equation which are given in terms of the CKVs of the space instead of restricting ourselves 
to the Noether point symmetries only, as has been done in fl3j . 

Below, we study the Lie point symmetries and the WdW equation for the potentials of the form (fl3l) which generalize 
the work done in si. 


4. LIE POINT SYMMETRIES OF THE WHEELER-DEWITT EQUATION 

We are considering a scalar field cosmological model which contains a quintessence scalar field with the potential of 
Ea. (fl3l) and a perfect fluid with equation of state parameter w rn = (7 — 1). Under these assumptions the Lagrangian 
of the field equations ED becomes 

^a, a, <p, <pj = —3 aa 2 + — a 3 </> 2 — Vb « 3 [a cosh (jxp) + /3 sinh (pcp)] q — Pmoa -3 ^ 7-1 - 1 . 


L 


(32) 
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From previous work on Noether point symmetries in scalar-field cosmology fl9| . and on dynamical symmetries, jl9l ] 
we know that this Lagrangian admits conservation laws when: a. the potential reduces to the exponential potential 
i.e. (3 = ±cc, and b. we have the so called Unified Dark Matter (UDM) potential (see Paliathanasis et al. 0 and 
references therein), i.e. p = ^ and q = 2, when the extra fluid is dust, namely ( w m , 7 ) = ( 0 , 1 ). 

In the following we consider q^/ 3 which implies that the current analysis generalizes the previous works of |3j,10 ■ 
The Hamiltonian (fTSl) of the field equations for the Lagrangian (1321) in terms of the momenta p a , and p 4 , is 

- + ^3 Pi + (^oa 2 3 [a cosh (pcj>) + fi sinh {pcj))] q + p m 0 a _3(7_1) ) = 0, (33) 

and the WdW equation m is 


- l^a^ ,aa + 9^3 ~ 12^2' I, ’“ + (^ 0 ° 3 cosh (W>) + /2 sinh {p4 >)] 9 + Pmoa 3(7 1} ) ^=0. (34) 

Applying the results of [0 , we find that the second-order partial differential equation (15H) admits the generic Lie 
point symmetry vector 

X = aX 1 +/3X 2 + a 0 ^d^ (35) 

where ag is a constant 2 , 


X\ = a 


3p, 


X2 — Cl 2 


~^r a sinh d a + cosh 84, 

—^r a cosh | d a + sinh | ^-P<f>j 84 


and the constants p , q, 7 are related as follows: 


4 

P = —rP- > q = - 2 , 7 = ^ + 2 . 

4 p 


That is, the effective potential of the field equations is 
Veff = I Vo a 3 


a cosh ( X-n(h ] + B sinh ( ^—pd> 

V 4 J V 4 > 


PmOCl 


—3(M+1) 


(36) 

(37) 

(38) 

(39) 


Therefore, for p = —1, we have that q = 2, 7 = 1 Le. we have the UDM potential with dust (for the exact solution 
and the observation constraints of that model see 0 ). 

If the perfect fluid is a barotropic fluid, that is the barotropic index 7 £ [1,2] then, from (1381) . p £ [—1,0) since 
p 0. However, if we require the perfect fluid to have a negative equation of state parameter, like a cosmological 
constant, then 7 £ [0, 2) which means that p £ [—2,0). Furthermore, when p m 0 = 0, i.e. there is no extra fluid, we 
have that p £ M.*. In the following, we apply the Lie symmetry vector (1351) in order to construct the invariant solution 
of the WdW equation (l34l) and to solve the null Hamilton-Jacobi equation of the Hamiltonian (l33l) in order to reduce 
the order of the field equations. In the following section we study the case a/3 = 0 and in appendix [B] we present the 
general solution for a/3 ^ 0. 


5. EXACT SOLUTIONS FOR THE cosh /sinh POTENTIAL 


In this section we determine the exact solution of the field equations and of the WdW equation for the quintessence 
scalar field. We consider the case a = 1, /3 = 0 (the case a = 0, f3 = 1 is equivalent to that case). Under the 
coordinate transformation 3 : 


a=(x 2 — y 2 ) 3l ‘ , (/) = — arctanh (-) , 

Ofl \ X / 


(40) 


2 The fields Xi,X 2 are CKVs of the minisuperspace (f20t . 

3 We assume fi < 0. However when p m o = 0 it is possible to have fi > 0. In that case all calculations remain valid provided we replace 

fi = —v in (1391) and in the subsequent coordinate transformations. 











the effective potential (f39|) becomes 


l + M 


V eff = V 0 (x 2 - y 2 ) » x » 2 


(41) 


and the WdW equation is 


(* 2 - y 2 ) 


2 \^ + l 


^,yy ^,xx 


(zVqX * 2 + 2p m ^j 


vh 


= 0 


(42) 


where Vq = |/i 2 Vo, p' m0 = | p 2 p m o • In these coordinates the Lie point symmetry vector (1551) is X = d y + a 0 ’P<9^. 
Therefore, the solution of equation (T751) admits an oscillatory term, i.e. (x, y) = e a ° v Q (x) where 


cf> - 

*' .cr.cr. 


(2V^x-^- 2 + 2p / mO + a 2 )^ = 0. 


Furthermore, in the WKB approximation, *F oc e lS and equation (l42l) becomes 


(x 2 -y 2 y +1 


dS 


dS 


-d-y) -{jTx) 


= 0 , 


which is the null Hamilton-Jacobi equation which describes the field equations. The solution of (1441) is 


S(x,y) = ay ± I \Jc\ + 2p' m o + 2H C 




(43) 


(44) 


(45) 


Therefore, the field equation is reduced to the following two-dimensional system, 



2VgX^^~ 2 


y 2 ) 


-(-+ 1 ) . 
u ) y = c 1 . 


(46) 


In order to simplify the system (TTBT) further we apply the transformation dr = (x 2 — y 2 ) ' 1+1 dt = a 3 ^ +1 ^dt, and 
the dynamical system becomes 

x' =T\Jc{+ Zp'mo + 2 VqX~^~ 2 , y' = a (47) 

The exact solution of the system (f47T) is expressed in terms of elliptic functions 4 . 

We perform a numerical integration of the non-linear system (14611 and in fig. [I] we give the evolution of the equation 
of state parameter for the scalar field w$ (a) and for the total fluid wtot (a) for various values of the constant Ci in 
the case p = —1. For p = —1 the extra perfect fluid is dust, i.e. (w m , 7 ) = (0,1). Concerning the values of the 
cosmological parameters, we use the Planck priors [50f, namely f2 m o = 0.308 and Hq = 67.8 km/s/Mpc which imply 
PmO = 3fl m o-£/o — 4.25 x 10 3 , in units of k = 8irG = 1. From figHJ we observe that the scalar field mimics the 
cosmological constant for small values of the constant ci, however for large values of C\ the scalar field has an EoS 
parameter > — 1. We find that within a physical range of the above cosmological parameters the corresponding 
dark energy EoS parameter deviates by ~ 1 — 2%. Furthermore, from the evolution of Wtot(ci) we see that there is a 
matter-dominated epoch. However, as the parameter ci increases, this epoch has shorter duration. In what follows 
we study the case ci = 0 and express analytically the scalar field and the Hubble function in terms of the scale factor. 


5.1. Subcase c 1 = 0. 


When ci = 0, from (147)1) we have y ( t ) = yo , hence from (14U1) it follows that x 2 = y 2 + a 3 ^. Furthermore, from the 
transformation (147)1) and from (147)1) we find that 


vl 

ifer + (% 2 + a- 3M ) 2 1 

— a 3ai (yo + a 3/i ) " 1 

vl 

j^f + (% 2 + a— 3/ T^ _1 

+ a~ 3 ^ (jjq + a~ 3fl ) M 1 


( 48 ) 


It is easy to see that when /x = — 1 then dr = dt, which is the UDM solution for UJ 2 = 0 of [19 
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w^(a) evolution forp= -1 (Dust) w (ot (a) evolution foru= -1 (Dust) 




FIG. 1: The evolution of the equation of state parameters w$ (a) for the scalar field and Wtot (a) for the total fluid for the 
Lagrangian (1321) with a = 1,(5 = 0. In (1461) we have taken the minus sign and p = —1 (dust fluid: Wm = 0). For the numerical 
solution we use (®,j/) t _ >0 = (0.01, 0.0099) , log(Vo) = 4, p m o = 3Q, m oHo with fi m o = 0.308 and Ho = 68Km/s/Mpc (or 
pmo = 4.25 x 10 3 in units of 87rG = 1). The solid line is for log(ci) = 0; the dashed line is for log(ci) = 1; the dotted line is for 
log(ci) = 2. 


where we have set 


0 _ Vo n — P m0 

“AO — o t T 2 > ^ L rnO — 


■m 


(49) 


Here, we would like to note that we call Hao the density parameter of the cosmological constant-like term and Q m o 
the parameter of the perfect fluid 5 . 

Therefore, for the scalar field density p$ (a) we have: 


P<t> (^) — ^QaqHqCl 


2-6 


Vo 


f^AO 


+ (y 2 + a~ 3 ^) " 1 ) + a _3/i (j/g + a~ 3fx ) 


-4-1 


Then © implies that 
r 2 


E 2 (a ) = = H m0 «- 3 ^+ 2 ) + H A 0 a - 6 (y 2 

-“o V 

We note that if yfi + a~ 3fl ~ a -3 ^ then 


^mO 


^A 0 


+ (Vo + «" 3m ) 


-T --1 


(y 2 o 


E 2 (a) — fl m oa + fl A0 ( 1 + yg 


^ a" 6 + a 3M 
^A 0 


(50) 


(51) 


(52) 


If 3 /o = 0, then (1521) becomes E 2 (a) = fl m0 a~ 3 ^ +2 ^ + fi A o which is obvious because when y 0 = 0, we have <j> = 0 
and V (</>) = Vo, which means that the scalar field acts as a cosmological constant. Furthermore, from (l5ll) and 
H (a = 1) = i?o, we have the constraint 


(l + 2/o) + f^AO (l + ( 


-1 = 0 . 


(53) 


In the following section we consider special values of the barotropic constant 7 = p + 2. 


5 In general, for yo 7 ^ 0? hold V2 m o + ^AO 7 ^ lj however the equality holds only when the constant yo = 0 = Vo] which means that 

the scalar field act as a cosmological constant, i.e. = —1, see eq. (153 D . 
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5.1.1. Dust fluid versus the effective dark energy EoS 


When the perfect fluid is dust then fj, = —1, 7 = 1 (w m = 0 for other cases see appendix A) Ea. fldTl) takes the 
following form 


E (a) — fi m O a + ^AO 


1 + 2 j/o a 3 


■Vo 


V ^AO 


+ Uo ) a 6 — 3 + A H(a) . 


(54) 


It should be mentioned that the last term AH (a) of the normalized Hubble function (1M1) introduces a cosmological 
constant-like fluid, dust and stiff matter. Furthermore, from (15311 we have the following algebraic equation 

n AO yo + (2Hao + £2m 0 ) Vo + (^m0 + ^ao — 1) = 0 (55) 

hence, the discriminant of the polynomial (1551) (for y%) is 

D = (2Hao + f2m 0 ) + 4 Hao (1 — H m o ~ Hao) > (56) 

and D > 0 when (1 — Q rn o — Hao) > 0. Recall that O m o € [0,1], Hao G (0,1], and because y$ > 0 we have the solution 

\j (Amo) 2 + 4 Hao — (2Hao + ^mo) 


vl = 


2 Hao 


(57) 


Let us now compute the effective dark energy EoS w^^ff for the scalar field model introduced above. It is well 
known that one can express the effective dark energy EoS parameter in terms of the normalized Hubble parameter 

El 


_i _ 2 d In E 

3 d In a 


_ 3d Inc 

WM ’ {a) - l-S!„,(o) 


(58) 


where fl m (a) = ■ Inserting the second equality of Ea. (l54l) into Ea. (l58l) . the effective dark energy EoS 

parameter takes the following form (see 0 ): 


w M f(a) = - 1 - 


1 din AH 
3 d, In a 


(59) 


which implies that any modifications to the effective EoS parameter are included in the second term of Eq. (El. 
Inserting Eq. (1+fl) into Eq. (15U1) it is straightforward to obtain a simple analytical expression for the effective dark 
energy EoS parameter: 


Wf>,eff(a) = -1 + 


2yla-* + 2yl^ + yl)a 

1 + 2 yla-z + yl^+yl) 


2 ia- 6 ' 


(60) 


5.1.2. The total case 


Our dynamical system is integrable in the case of a single perfect fluid. Here, we introduce two perfect fluids 
(for example dust and radiation). Consider that we have dust (w m = 0) and another perfect fluid with equation of 

state parameter Pf = Wfpf (for radiation Wf = 1/3), and ("C 1. The latter implies that the equation of state 


parameter which is associated with the two perfect fluids is 

_ _P m + Pf_ WfPf _ w f{ nt) 


Pm+ Pf Pm+ Pf 1 _|_ 


Wf 


12 / 
1 I'm 




(61) 


that is, 7 = 1 + w m and p = — 1 + w m . When Wf > 0 we have that w m > 0 and when wf < 0 holds we have w m < 0. 
We replace m in m and perform a Taylor expansion near w m = 0(7 =lor/x = —1). We find 


E 2 (a) = -E 2 =1 (a) - 3H m0 a 3 In (a) w m + VL ao F (a) w m 


(62) 
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w (a) evolution for u=-1 ±0.05 

tot' ^ 



FIG. 2: The evolution of the equation of state parameters Wtot (a) for the total fluid for the Lagrangian (1321) with a = 1, fi = 0 
where in ([46} we take the minus sign and ci = 0. For the numerical solution we use (x , y) t _ >0 = (0.01, 0.0099) , log(Vo) = 4 
and pmo = 4.25 x 10 3 . The solid line is for p = —1 [dust: (ui m , 7 ) = (0,1)], the dash-dotted line is for p = —0.95 and the dotted 
line is for p = —1.05. 


where the normalized Hubble parameter E^ = 1 (a) is given by Eq. (l54l) and 

F (a) = 2 In + a 3 ) (l + y% a -3 ) 2 - 6 (j/q + a 3 ) a ~ 3 In a 

where F (a —> 1) = 2 In (j/q + l) (l + y q ) 2 and when y 0 = 0, F (a) = 0. In fig. © we give the numerical solutions of 
the total EoS parameter wtot (a) for p = 1 ± 0.05 and ci = 0. 

In appendix El we give the exact solutions for the Hubble function, FI (a), for the cases where 7 = | (radiation 
fluid) and 7 = | (curvature-like fluid). 


6. DYNAMICAL ANALYSIS 


In order to complete our analysis of the model with Lagrangian (1321) . we perform a dynamical analysis of the field 
equations by studying the fixed points of the field equations. We introduce the new dimensionless variables @, [45[ 


W p 

y = ^ • Urn = 


— \ = -Y± 

’ " m 3H 2 V 


VQH’ y V3H 

and the lapse time N = In a. In the new variables the field equations reduce to the following first-order ODEs 

^ = -3x + + Y [(1 - Wm) X 2 + (1 + Wm) (l ~ 2/ 2 )] 


(63) 


(64) 


= ~Y- Xx y + \y [(f - W ™) X 2 + (1 + W m ) (1 - y 2 )] 


dX 


— = - V 6 A 2 (T-l)x 


(65) 

( 66 ) 


where T = and the Friedmann equation ([ 6 ]) gives the constraint = 1 — where 


v^v 

Vl 
>0 

In this case the second Friedmann equation 0 becomes 

2 H 


= x 2 + y 2 . 


3 H 2 


= -1 - W m - (1 - Wm) X 2 - (1 + Wm) (l ~ y 2 ) 


(67) 
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TABLE I: Fixed points and cosmological parameters 


Point(x,y, A) 


0 

(0,0, A) 

A ( ±) 

(1,0 ,±qp) 

B (±) 

(-1,0, ±qp) 

C 

(0,1,0) 

D (±) 




W tot 


Acceleration 

1 

Wm 

$ 

Wm 3 

0 

1 

1 

No 

0 

1 

1 

No 

0 

-1 

-1 

Yes 

0 

1 1 ( qp ) 2 

1 ^ 3 

1 1 (qp ) 2 
' 3 

\qp\ < a/2 


E (±+) 

( 1 \/6(l+'ui m 

^ 2qp 


, 3(l+™ m ) 

(qp)* 

Wm 

E (±~) 

( 1 A/6(l+UJ m ) 

^ 2 qp 

,-^pE, ± „ p ) 

1 3(1 + ^771.) 

(qp) 2 

Wm 


Wm Wm O g 

Wm Wm O g 


TABLE II: Eigenvalues of fixed points 


Point 

mi 

m2 

m 3 

O 

0 

| (1 + Wm) 

-| (1 - W m ) 

A ( ±) 

3 (1 - Wm) 

±2pVE 

3 T#9P 

B (±) 

3 (1 - Wm) 

=p2pv/6 

3 ±^q P 

C 

-3 (1 + Wm) 

“I (l- V 1 ” 4 ® 2 ) 

~f (t + V 1 - 4 ip 2 ) 

D (±) 

-3 (1 + Wm) + (qp) 2 

-3+^|l! 

2 qp 2 

E (±,±) 

f (1 + W m) 

“I [OL-WrrO + ^l 

“I [(!-«» m )- £] 


which gives that the total EoS parameter wtot as a function of w m , x and y: 


Wtot = Wm + (l - Wm) X 2 - (1 + Wm) y 2 ■ (68) 

Furthermore, the EoS parameter for the scalar field is = x x 'i+ y ‘± . Note that at any point (xo,yo. A), from (167)) 

2 

the solution of the scalar factor is a power law as long as Wtot = const.; that is, a ( t ) oc 1 3 < 1 +™tot) for Wtot ^ —1 and 
a (f) = aoe Hot for w to t = -1 • 

In the following we consider in (ld2l) /3 = 0, so that the potential of the scalar field is V (</>) = Vo cosh 9 (p(j>) (46l - l4g| . 
For this potential we write T (</>) as a function of A, i.e. T (A) = 1 + ^and equation (1551) becomes 

^ (qp - X) {qp + X) x. (69) 

Equations (1641) . (1651) and (1691) describe an autonomous dynamical system in the E 3 space. Furthermore from the 
constraints 0 < < 1, y > 0, the variables (x, y) are bounded in the ranges x G [—1,1], y G [0,1] from which follows 

that that the points (x, y) belong to a half disk; however for the parameter A there is no constraint that implies that 
A G R [ill dll- Furthermore, we consider w m G (—1,1). The fixed points of the dynamical system (15T1) . (1551) and 
(1691) and the corresponding cosmological parameters are given in Table |U The eigenvalues of the linearized dynamical 
system near the fixed points are given in table [Tip. 

Point O exists for all values of the parameter A and corresponds to the matter epoch (fl m = 1); the total EoS 
parameter is wtot = w m ■ Since there exists at least one positive eigenvalue, m 2 > 0, the point O is always unstable. 
At this point the universe accelerates if and only if w m < — |. At the points A(±) and .B(±) the universe is dominated 
by the kinetic energy of the scalar field (fl m = 1, V ((/>) = 0) which means that the scalar field acts as a stiff fluid, 
i.e. ptj, oc a~ 6 which provides a decelerating universe. These points exist when A = ±qp, for arbitrary q,p. For these 


Where in table EDa 



Wm ) ^24 (1 + Wm ) 


(7 + 9 Wm) (gp) 2 ) 
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FIG. 3: Phase portrait in the x — y plane and in the E 3 space for the potential V (</>) = Vo cosh 9 (p(j >). Left-hand figures are for 
the variables ( p , q, Wm) = (1, —1.5, 0) and the right-hand figures are for the variables ( p , q, w m ) = (1, —3, 0) . For ( pq ) 2 > 3, the 
stable points are the points -E(±) (scaling solutions) while for (pq) 2 < 3 the two stable points are the points D±. The solid lines 
are for the initial condition A = pq, the dashed lines for A = — pq, and the dotted lines for A = 0. 


points there exist positive eigenvalues of the linearized system, mi > 0, for w m G (—1,1), hence these critical points 
are always unstable. 

Point C is the de Sitter solution, (f2 m = 0, w to t = — 1) where the scalar field acts as a cosmological constant and 
the matter component vanished. This point exists for all values of the constants q,p and could be the future attractor 
of the universe. From the eigenvalues of Table [TT] for that point we have that it is stable when q > 0 (for a similar 
solution see [HI). The points Dp n correspond to a scalar field dominated universe (O m = 0) and exist only when 

\qp\ < a/ 6- The total EoS parameter is that of the scalar field = — 1 + which gives an accelerated universe 
when \qp\ < y/2. The points -D(±) are stable for q < 0 and \qp\ < ^/3 (1 + w m )- Hence, we see that D(±) are stable 
points and describe an accelerated universe when | and \qp\ < \[2. Furthermore, in the limit q —> 0 _ , these 

points correspond to the de Sitter universe; when q —> 0“ then V (</>) —> Vq. 

Finally, the points i/±,±) are the so called ’scaling’ solutions where f l m = 1 — and the scalar field mimics 

the matter component of the universe, i.e. vu$ = w m . The points E/± i+ \ exist when qp > -^/3 (1 + w m ) and they are 
stable when q,p < 0 whereas the points E(± exist when qp < —1/3 (1 + w m ) and are stable when q < 0, p > 0. 
The total EoS parameter is Wtot = uim so they lead to an accelerated universe when w m < — -. 

In Eg. [3j we give the two-dimensional phase portrait in the x — y plane and the three-dimensional phase portrait 
of the model with values ( p,q,w m ) = (1,—1.5,0) and (p,q,w m ) = (1,—3,0). We observe that for (pq) 2 < 3 the two 
stable points are the points -D(±) whereas for (pq) 2 > 3 the stable points are Er± t _y 

It is important to study the case when the constants q,p,w m are related to the constant fi so as to render the 
integrable field equations (l38l) . This is the case we studied in the previous section. Since we considered w m £ (—1,1) 
we have that fi £ (—2, 0). 

Hence, for the integrable case we have that the points O , B/±\, D( ± ) exist and they are always unstable. 

The point O has wtot < —| when p G (—2, — |) and the points -D(±) describe an accelerated universe so long as 
p £ ^—2, —2 + ■ The point C exists and it is the unique stable point. Finally, the tracker solutions, i.e. points 

E(±,±), do not exist for p G (—2, 0). The existence and the stability of the fixed points for general values q,p and for 
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TABLE III: Fixed points and their stability for the general potential and for the integrable subcases 


Point 

Existence 

Stability 

Stability for y £ (—2,0) 

Acceleration 

O 

p,q el* 

Unstable 

Unstable 

y£ (-2,-|) 

A (±) 

p,?el* 

Unstable 

Unstable 

No 

B{±) 

p,q 6 R* 

Unstable 

Unstable 

No 

C 

p, q £ R* 

Stable for q £ R* + 

Stable 

Yes 

T>(±) 

\qp\ < V6 

Stable for q £ R*~ , \qp\ < ^/3 (1 + w m ] 

Unstable 

y£ (—2,-2+2^2) 

E (±,+) 

qp> y/3 (1 + w m ) 

Stable for q £ R*“ , p < 0 

$ 

$ 

E (±,~) 

QP < — \/3 (1 + w m ) 

Stable for q £ R*“ , p > 0 

$ 

$ 


Phase portrait in the x-y plane for p=-1 (Dust) 

Phase portrait for 

=-1 (Dust) 



c 


, - 

/ / \\ 
/ / 

/ / / Vv 

/ / vlV 

\ 

AV X \ 

\ 

/ / » 

^\/r\ \ \ 

/'' / 

ivA \ \ V 

/ / ' 
l / , J 

lift'. \ V \ \ 

\ X \ \ 

fl lv \\ \. \ \ 1 

/|'W\ Nr. X \ 

'/V v \ 



FIG. 4: Phase portrait in the x — y plane and in the E 3 space for the potential V ( <j >) = Vb cosh 9 (pep) with variables ( p , q, w m ) = 
^--^,2,0^ , which corresponds to the integrable case for y = —1. The point C is the unique stable point. The points -D(±) 

act as attractors in the plane (x — y) for A = ±-^2. From the right-hand plot we observe that the solutions reach a boundary 
where fl m o = 0, > —1 and from there they move to the de Sitter point C. The solid lines are for initial condition A = pq, 

the dashed lines for A = — pq, and the dotted lines for A = 0. 


the integrable case are given in Table IIIII 

In fig. [3] we give the two-dimensional phase portrait in the x — y plane and the three-dimensional phase portrait of 
the the model with values ( p,q,w m ) = ^,2,0^ , which correspond to the integrable case for y = —1. The point 

C is the unique stable point. We observe that the points Z? (± ) act as attractors in the plane (x — y) for A = 
and the solutions reach the boundary where fl m = 0, and move to the de Sitter points (w^ = — 1). It is important to 
note that the existence of conservation laws in the field equations which follow from the Lie point symmetries of the 
WdW equation, i.e. the dynamical system is Liouville integrable, gives us constraints on the free parameters of the 
model so that there exists a unique stable point which describes the de Sitter universe. 


7. CONCLUSIONS 

We have applied Lie symmetry methods in order to extend the works of Rubano and Barrow fsf] and Paliathanasis 
et al. 0 in scalar field cosmology for a general family of potentials, V(4>). We have shown that there exists a unique 
connection between the Lie point symmetries of the WdW equation and the conservation laws of the field equations. 
We considered a general form of V(cf>) which contains hyperbolic functions for a scalar field with a perfect fluid and 
we have investigated the existence of Lie point symmetries of the WdW equation. This approach is more general than 
the application of Noether point symmetries. We recovered the result of Q, that is, that in scalar field cosmology 
amongst the variety of V (</>) potentials the hyperbolic types play a key role because they admit conservation laws. 
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Moreover, based on the Lie point symmetries of the WdW equation, we have obtained the exact solutions of the field 
equations. Finally, we have performed a dynamical system analysis by studying the fixed points of the field equations 
in dimensionless variables. We found various dynamical cases among which, if the current cosmological model is 
Liouville integrable there is a unique stable point which describes the de-Sitter universe as a late-time attractor for 
the dynamics. 
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Appendix A: Solutions with radiation and curvature 


Here we provide some details concerning the solutions of section [All Specifically, for a radiation perfect fluid, 7 = 
(hence p = — | (w m = 1/3) and f2 r0 = fi m o) Ea. rtbTl) gives 


E 2 (a) = 


H 2 (a) 

H 2 0 


— Q r oci 4 + Hao 


1 + 3 yla 2 + 3 y„a 4 + yl 


4-4 


( ^r-0 2 \ -6 

n7 + !to *“ 


(Al) 


On the other hand, when 7 = | (or \x = — |), the perfect fluid has an equation of state w m = — ^ (which can also 
be seen as the curvature term in a non-flat FLRW spacetime); then fix 0 = fl m o and: 


E 2 (a) = 


H 2 (a) 

Hi 


= fix 0® 2 + f^AO 


Vo + a 4 a 2 + yl 


( 


ft 


K 0 


V ^AO 

This is a solution of the scalar field cosmology in a curved FLRW spacetime. 


+ \JVo+ a « 


— 6 


(A2) 


Appendix B: Exact solution for a general potential 

In the general case where a ^ 0 and /3 = 1 we apply the following coordinate transformations, 


a = 




y 


2 V 6 
3 p 


arctan h 


y 


In the new coordinate system the WdW equation becomes 

1 + U- 


HI) 


y 


-ll-i 14- 

a z 


,xx + —^,xy + 4-,yy + e 2 x V 2 + 2p' m ^j T 


= 0, 


(Bl) 


(B2) 


where V/ = |^ 2 a 5* 2 Vo, p' m0 = §/r 2 p m o and the Lie symmetry vector is (1551) X = d y + a^dy. Therefore, the 
invariant solution of the WdW equation (IB2I) is -F (x,y) = e Qoy <E> ( x ) where ( x ) satisfies the following second-order 
ODE 


- 1 - 


1 


$ 


++ (2Vq x - 1 2 + 2p(„o + ao) 


$ 


= 0 . 


When a = 1, which is the case of the exponential scalar field, from (IB3|) we have 


^, a | (^PmO + a o) 

<P (a:) = exp--- -x ■ 

2 ao 




a ° (l + X ) 


-±-l 
X * 


(B3) 


(B4) 


Furthermore, in the WKB approximation the WdW equation (IB2I) becomes the null Hamilton-Jacobi equation 




l+M r 


- 1- 


1 \ fdS_ 
dx 


r,\ 2 


2 /as\ /as 
a \dx J \dy 


+ ( ) + ( 2V ° ,;C * 2 + 2f> ' m0 ) 


= 0 (B5) 
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with solution of the form S (x, y ) = Si (x) + C\y, where 

Si (x) = =F — a 2 J ( Cl + \! 2V o (« 2 - !) x~^~ 2 + a 2 cf + 2 (a 2 - 1) p' m ^j , \a\ ^ 1 (B6) 


and 


Si (x) = 


-YL x -$~* - i Cl - 

ci 2 ci 


dx , a = 1. 


(B7) 


From the solution of the Hamilton-Jacobi equation (1B5I) we can reduce the equivalent Hamiltonian system of the 
field equation to the following system of first order equations 


x = 


y = 


- i- 


l 


~Px + Py 
a 


1 

Px + ~Py 
a 


(-!) 






l + M 


(B8) 

(B9) 


where p x = and p y = . 


dx t'v ~ d y - 

In order to make the reduced system simpler, we apply the transformation dr = a _3 (/ i + 1 ) dt, and the system 
becomes 


x' = - 1 - 


/ 1 \ 1 ,1 

1 - 2 )P* + ~Py > y = “P® + Py • 

\ cr / a a 


(BIO) 


For the exponential potential (a = 1), from system (1B10I) we find the solution in closed form. The solution is 

x(t) = cit + cq, (Bll) 


y( T ) = 


Vo 




(ciT + Co) " 1 - —T + ^-r + 2/0, 

Ci 2 


(B12) 


where the scale factor is a (r) = (x 2 (r) — 2x (r) y (t)) 3,i . This is the solution of the exponential scalar field with 
matter in the Einstein frame. Recall that the matter has an equation of state parameter of the form w m = fi + 1. 


1. Special solution with p m o = 0, ci = 0 


When p m o = 0 and Ci = 0, from the dynamical system (IB10I) we have the solution 


■(t) = 




+ 1 X 0 T + X 1 


2(l + /i) 


, y( T ) = — 2 — ~y°- 

or — 1 


where x$ = ^^/2 Vq (a 2 — 1), and e = ±1. 

In the case of /i = — 1 the solution of the system (IB 101) is 


x (r) =x 1 e X0T , y (t) = — 


a 2 — 1 


x (t) - y 0 


(recall that for p = —1, we have dt = dr). Hence, the solution for the scale factor is 

a 2 - 1 


(t) = 


Q^l „2x n t 


a 2 — 1 


a 


2 Vo 


(B13) 


(Bll) 


(B15) 


Furthermore, from the singularity condition a (0) = 0, we have that yl = which implies a (t) = 


a 1 (e 2xot — l) 3 , where oq = We obtain H (t) = ^ jT 0 t\ , and f (a) = ^ln 1+ 

fore, we can express the Hubble function in terms of the scale factor, i.e. 


E 2 {a) = 


H 2 (a) 

Hi 


— ^A0 T “1“ ^s/0^ 


-6 


. There- 


(Bi6) 
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where 


O 4x ° 

fiA0= 9i7f’ 


n 


- 8 x o n 3 

m0 ~ 9 Hf 1 


Q 


sfO — 


a 6 

9tf 0 2 1 


which means that the scalar field introduces an effective dark matter component, namely Q m (a) = Q m a 3 /E 2 (a) in 
the cosmic expansion. 
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